Instanton interaction in de Sitter spacetime 



Dimitrios Metaxas 



Department of Physics, 
National Technical University of Athens, 
Zografou Campus, 15780 Athens, Greece 

metaxas ©central, ntua.gr 



Abstract 

Because of the presence of a cosmological horizon the dilute instanton gas 
approximation used for the derivation of the Coleman-De Luccia tunneling 
rate in de Sitter spacetime receives additional contributions due to the finite 
instanton separation. Here I calculate the first corrections to the vacuum 
decay rate that arise from this effect and depend on the parameters of the 
theory and the cosmological constant of the background spacetime. 



1 Introduction 



The calculation of the decay rate of metastable vacua in flat and curved 
space-time [1, 2] and zero or finite temperature [3] is of great interest due 
to various applications to problems at the interface of particle physics and 
cosmology. Especially, the recent discovery of a multitude of string vacua 
suggests that a detailed understanding of these rates for different parameters 
of the models may be relevant to the investigation of the various cosmological 
phase transitions that took place in the evolution of the universe as we know 
it [4]. In particular, it is important to have quantitative results for the 
dependence of the decay rate on physical parameters such as the temperature 
[5] and the Hubble expansion rate or the curvature of the background space- 
time [6]. 

To this effect, the usual semiclassical calculation at zero temperature in 
flat space-time, which I review in Sec. 2, needs various modifications and 
admits additional contributions depending on the physical situation at hand. 
In Sec. 3 I proceed to examine the corrections that arise in de Sitter space- 
time because of the presence of a cosmological horizon. The dilute instanton 
gas approximation then receives corrections, the first of which arises from the 
two-instanton interaction, and I calculate it in the thin wall and fixed back- 
ground approximation. The fact that the instanton interaction in de Sitter 
space-time turns out to be repulsive, essentially because of the cosmologi- 
cal expansion, and unlike most instanton and soliton interactions in scalar 
field theories, makes the calculation feasible, without the need of additional 
cutoffs. 

The corrections obtained here are, of course, subleading with respect to 
the semiclassical exponential factor, they are of interest, however, because 
of their dependence on the various parameters of the theory, especially the 
cosmological constant of the background space-time. Finally, in Sec. 4 I 
conclude with some comments regarding possible extensions of this work. 

2 Review of the flat space-time results 

In order to calculate the false vacuum decay rate for a scalar field theory in 
flat, 3 + 1-dimensional space-time and zero temperature one evaluates the 
imaginary part of the energy of the false vacuum, obtained from the path 
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integral 

Z(T) = y[#]exp(-5 B (0)), (1) 
with the Euclidean action 

r r T / 2 

S E ((f>) = d 3 x dt E L E , (2) 

J J — T*/2 



b = ^(^) 2 + W) (3) 



and 

1 

2 

where the potential U(<f>) has a metastable minimum at <fi = <pf and an 
absolute minimum at = <f> t , and the path integral is restricted to paths 
satisfying 0(-T/2) = <j>(T/2) = <p f . 
Then, in the limit T — > oo, 

Z{T) w exp(- J B / TO), (4) 

where f2 is the volume of space and the energy density of the false vacuum 
and the imaginary part of this quantity will give the decay rate 

r = -2 ImlnZ ' r », (5) 

TO. y ' 

that is, in this case, the bubble nucleation rate, or the number of true vacuum 
bubbles produced per unit volume per unit time. 

The path integral can be evaluated in the saddle-point approximation by 
obtaining an instanton (bounce) solution, <p b , of the Euclidean field equation 

0*= w (6) 

Then 

Z(T) « Z{<j> f ) + Z{4> h ) + Z 2 + ... (7) 

where the first contribution comes from a homogeneous false vacuum config- 
uration, 

Z{<t> f ) = [det S£(0/)]- 1/2 exp(-5 E (0 / )) (8) 

(S'e denotes the second variation of the Euclidean action), the second con- 
tribution to (7) comes from the bounce configuration, 

Z{<t> b ) = '-QTJ | det'S'^fo) r 1/2 exp(-S E {<k)) (9) 



and higher contributions come from multi-bounce configurations. In the 
evaluation of (9) one takes into account that around the bounce configuration 
S'e has one negative mode and four zero modes (which are excluded from the 
determinant). The treatment of the negative mode gives the factor of i/2, 
the zero modes can be treated by changing to collective coordinates, which 
gives a Jacobean J, and the integration over these coordinates gives a factor 
of QT, the total physical volume of the system. 

The integration over the zero modes can be better illustrated in the case 
of the + 1-dimensional case of a simple quantum-mechanical system with 
the Euclidean action 

SWtf) = dt E Q(^) 2 + [/(</>)) , (10) 

with a similar potential that admits a (one-dimensional) bounce solution, 
(f>b~\t), to the Euclidean equation <9 2 = dU/d<f). Because of translational 
invariance there is a single zero mode. 

4> a = $\t-8) (11) 

is a one-parameter set of solutions with the same action, and O = d t 4> s / \J 
where is the action for <$\ is the corresponding zero mode. 

The corresponding path integral can be evaluated in the semiclassical 
approximation around the saddle point [7], after writing cf>(t) = <p s + ip, 
expanding ip = Y, a n^Pn, in eigenfunctions of Q — 5 2 S ( - lS) /5 2 4>, such that 
Qip n = \ n ip n , and inserting the identity 

1 = A(0) J ds6[<<l> , </>-</>,>], (12) 

where < f,g >= J f(t)g(t)dt. The action to second order is S — + 
| X) AnQ; 2 , and the divergence of the integration over the zero mode is taken 
care of by the delta function, resulting in an overall factor of T, the "volume" 
of the system, and a Jacobean factor, J = A(<j>) A(0 s ) in the semiclassical 
approximation. 

Going back to the 3 + 1-dimensional case, one continues with the higher 
contributions to (7), which come from multi-instanton configurations: for 
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example, for the two-bounce configuration, dividing the total space-time vol- 
ume into two volumes V\ and V2, each containing one instanton, one obtains 
[8] 

Z 2 = |[dfl e -*-*- = ±J[d<p\ Vl j[d<P\ V2 e-k LE e-^ 2 L * (13) 

1 Jld^ve^ J[#],e-* 



2!/[d0] V2 e-So/[d0] Vl e-Sb 
/ ^ \ 2 



7, 
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and 



= */ + */g + *4(f*)' + " (M) 

From (8, 9) we get 

I* = \-VtTAe~ B (15) 

where A denotes the Jacobean and determinant factors and 

B = S(<j> b )-S(<f> f ), (16) 
and, finally, from (5) we have 

T = Ae- B . (17) 

As a final note one should mention that the bounce solution of (6) is ob- 
tained in flat space-time using the fact that it enjoys 0(4) symmetry [9], and 
that the dilute instanton gas approximation that was used when considering 
widely separated instantons in (14) can be justified when B » 1. 

3 The bounce and the instanton gas in de Sit- 
ter space-time 

When extending the previous results to include the gravitational effects on 
the false vacuum decay rate one assumes again 0(4) symmetry (which has 
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not been rigorously proven in this case) with the metric 

ds 2 = dr 2 + p(T) 2 dn 2 3 , 
and then the Euclidean action becomes 

3 



S fi = 2vr 



5* + 



8nG 



(pP + PP ~ P) 



for which the Euclidean field equations are: 



P 2 -1 



8ttG 



-P 



1 • 



<t> 2 - u{4>) 



( l ) +— ( t ) = -77 ■ 
p d<p 

I will work in the fixed background approximation where 

u(<j>) = u + u{4>) 



(18) 
(19) 

(20) 
(21) 

(22) 



with | U(4>) \« U , and a mass scale p such that ~ /i, U{<pf) ~ 

U{(f> t ) ~ Uo » p A . As before, 4>t and 0/ are, respectively, the absolute 
minimum (true vacuum) and the relative minimum (metastable vacuum) of 
the potential. 

Then one has a Euclidean de Sitter background spacetime which is topo- 
logically a four-sphere with 



where 



H 



8ttG, 



(23) 



(24) 



and the coordinate r extends from 0, on the "south pole", to jj, on the 
"north pole" of the sphere. 

The bounce solution to (21), in the thin wall approximation of [2], takes 
values close to <j> t at the "south pole", r ps 0, and falls rapidly to 0/ at larger 
values of r, satisfying the boundary conditions <fi = at r = and r = jf. 
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Provided that the bounce radius, Rb, is much smaller than one can 
justify the dilute instanton gas approximation described in the previous sec- 
tion and proceed to derive the Coleman-De Luccia bubble nucleation rate, 
?cdl = Ae~ B , where 



and A is the functional determinant pre-factor around the single bounce 
solution in de Sitter spacetime. Since, however, the Euclidean de Sitter four- 
sphere has a finite extend, r extends up to the de Sitter horizon 



and there are always contributions from the interactions between instantons, 
multi-bounce configurations, as they are distributed on the four-sphere, es- 
pecially when their radius grows larger. 

Here I will estimate the first corrections to Tqdl that appear because 
of this effect. In order to get a quantitative estimate I will assume that we 
are still inside the limits of the previous approximations (in particular, the 
condition B » 1 also has to hold, as in the flat case, and the thin wall and 
fixed background limits are also assumed) and obtain, therefore, a correction 
to the pre-exponential factor for T, which is of interest, however, since it 
depends on the parameters of the model and especially on the cosmological 
constant of the background de Sitter space-time. 

A first estimate of the corrections involved comes from considering the 
two-instanton interaction, as it appears in (14) and its generalization to the 
Euclidean de Sitter spacetime. The method involved is similar to the calcula- 
tion of soliton interactions in [10] and a basic reason that it can be used here in 
order to get quantitative results is that the instanton interaction in de Sitter 
spacetime turns out to be repulsive, essentially because of the cosmological 
expansion rate, and unlike most scalar soliton and instanton interactions in 
quantum field theory. The two-instanton configuration that gives the great- 
est contribution to the saddle-point evaluation of (14) will consist, therefore, 
of two bounce solutions symmetrically on the "north pole" and "south pole" 
of the de Sitter four-sphere. 

In order to justify the previous statements we consider, accordingly, a 
two-instanton configuration of the form 




(25) 




(26) 



(27) 
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where 

= - (f) f (28) 
and (j) b a solution to the bounce equation 

+ 3H cot (Hr)<j) = (29) 

When the two-instanton separation s is very close to the previous con- 
figuration is an approximate solution of the Euclidean field equations and 
one can estimate their interaction as a function of s. As in [10], the inter- 
action will come from the overlap of the "tails" of the two widely separated 
instantons, that is around the "equator" r = One needs, therefore, the 
asymptotic behavior of the bounce solution to (29). First we write around 
the false vacuum 

W) = f/ o + ^(0-0/) 2 (3°) 

with 

fi = ^. (3D 

After setting x = r — , writing cot(Hr) fa —Hx for small x, changing 

to = u(x)e^ x2 and neglecting terms of higher order in x in the final 
expressions, we find the asymptotic behavior for = — 0j around r ~ ^ : 



0^ C e"VM- (32) 

with c a constant of order \i. 

Now one can estimate the two-instanton interaction SB as the difference 
between the value for B of the configuration S and twice the value for a 
single bounce, 

B((j) s ) = 2B + 5B, (33) 
and, using the previous expressions, we find 



5B = 2tt 



fdrf 



^t)Us-t) + PH 2 Mt)Ms-t) 



(34) 



= (27rV / dr—^m\Hr)]-H 2 e- H ^ Z i s 

V Jeq H ) 2 



27r 2 c 2 ^--H 2 e- H V^ s . 
3H 4 2 
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Here, the integration inside the parenthesis in the second line of the previ- 
ous expression is done around the "equator" r « since, however, the 
integrand falls to zero away from this value, I have extended the region of 
integration to the entire range, the difference being a numerical factor of 
order unity that can, in any case, be absorbed in c. 

One can see immediately from this expression that the effective instanton 
interaction in the background de Sitter spacetime is repulsive, and, conse- 
quently, that the path integral is dominated by growing values of s ~ ^. 

If this instanton interaction term did not exist, one would proceed from 
expressions like (14) to derive Ycdl as before, the integration over s, from 
s = to s = jj, corresponding to a collective coordinate contributing to 
the total volume factor in (15). Now, however, this integration will give a 
different factor, leading to an additional contribution to the pre-exponential 
factor for the false vacuum decay rate. In the two-instanton interaction, the 
factor Z 2 in (14), after treating (27) in the manner of (12), gets modified to 

z *4 QZ '(t) 2 ' (35) 

with 

Jo dse 

a = n/H (36) 

and in the general n-instanton contribution to (14) one may consider the 
nearest neighbor interaction and integrate out each instanton to get the gen- 
eral factor 

*-4© Vl - < 37) 

Overall this leads to a correction to the pre-exponential factor of the Coleman- 
De Luccia result 

T = aYcoL (38) 

with 

a = - [ n dxexp(- ai e- a2X ) (39) 

7T JO 



where a\ = A ^f and a 2 = y(3 — |. When — << 1 this has the asymptotic 



expansion a ~ 1 — — , where 7 = 0.57721... is the Euler-Mascheroni 

constant. 
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4 Comments 



In summary, using the thin-wall, dilute instanton gas approximations of [2] 
and the fixed-background approximation as described before, I derived the 
first corrections to the pre-exponential factor of the Coleman-De Luccia vac- 
uum decay rate that originate from the instanton interactions inside the finite 
horizon of the background de Sitter spacetime. 

It should be stressed again that the general configuration considered in 
(27) is only an approximate solution to the Euclidean field equations, when 
the instanton separation is very close to the horizon value s = ir/H; because, 
however, of the repulsive nature of the instanton interaction, the integrals 
considered are dominated by these values (this is reflected by the fact that 
the leading term of the above expressions for the corrections is unity). I 
should also mention that the analysis presented here concerns essentially 
what would be one of the several zero modes, not considering the rest or 
the important question of the negative mode(s) [11]. Possible additional 
contributions arising from the treatment of these modes will be considered 
in future investigations. 

It would also be interesting to see whether some of the approximations 
used can be relaxed; it is possible that this may be true for the thin wall 
approximation, since the instanton interaction comes mainly from the overlap 
of the "tails" of two bounce configurations. It would also be useful, from 
the phenomenological point of view, to investigate the limits of the fixed- 
background approximation. Hopefully these issues, as well as the application 
of these results to models of the cosmological landscape will be the subject 
of future work. 
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